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Employing the operator algebra of the conformal group and the conformal Ward identities, 
we derive the constraints for the anomalies of dilatation and special conformal transformations 
of the local twist-2 operators in Quantum Chromodynamics. We calculate these anomalies in 
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the leading order of perturbation theory in the minimal subtraction scheme. From the conformal 
consistency relation we derive then the off-diagonal part of the anomalous dimension matrix of the 
conformally covariant operators in the two-loop approximation of the coupling constant in terms 
of these quantities. We deduce corresponding off-diagonal parts of the Efremov-Radyushkin- 
Brodsky-Lepage kernels responsible for the evolution of the exclusive distribution amplitudes and 
non-forward parton distributions in the next-to-leading order in the flavour singlet channel for 
the chiral-even parity-odd and -even sectors as well as for the chiral-odd one. We also give the 
analytical solution of the corresponding evolution equations exploiting the conformal partial wave 
expansion. 
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1 Introduction. 



Quantum Chromodynamics (QCD) — a relativistic local quantum field theory based on the colour 
gauge group S77 c (3) — is well established nowadays as a microscopic theory of strong interaction. 
However, due to presently unsolved bound state problem, which remains a challenge for future 
studies, it is not possible to predict measurable observables from first principles. The standard 
approach which overcomes this difficulty is based on factorization theorems |l] which allow to 
separate soft and hard physics in a given process. The soft part is described by expectation values 
of appropriate operators sandwiched between hadronic statesf] supplied with appropriate evolution 
equations which govern their scale dependence. Due to lack of any reliable non-perturbative 
machinery, the former part can be fixed only from experimental data. At the same time the hard 
subprocess and the evolution of the soft function can be systematically calculated with the help of 
perturbation theory. The success of perturbative QCD in predicting the momentum dependence 
of hadronic observables serves as the main and the most important argument for the correctness 
of the theory. Therefore, one is interested in increasing the accuracy of the theoretical predictions 
for their scale violation. 

In general the calculations carried out beyond leading order in the coupling constant are very 
difficult and require substantial computer power. There are continuous advances in higher loop 
calculations of inclusive processes. For instance, the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution kernels |2|, ||, §] for the twist-2 parton distribution functions (see Refs. || |] 
for transversity) are known in next-to-leading order (NLO) J7|, || [| and the calculations for the 
next-to-next-to-leading order (NNLO) are in progress. For exclusive processes no comparable 
technical effort was invested so far. 

For exclusive reactions at large momentum transfer the soft physics is contained in distribution 



amplitudes (DA) |L0], |Tl|, [12|] or more general non- forward parton distributions [13], |Tj, [15], [Tql , 
which depend on the momentum fraction x carried by the struck parton, and possibly other 
kinematical variables such as the skewedness parameter of the process. The scale dependence 
of the mesonic DA is governed by the Efremov-Radyushkin-Brodsky-Lepage (ER-BL) evolution 



equation |ll], |17|, [18 



d r 1 

f^-rt^fav) = / d v v ( x ,y\ a s{v))<f>{y,fj'), (i) 



dfi 2 ' Jo 

where the evolution kernel V(x, y\a s ) = ^V^(x, y) + (j^-) V^(x,y) + . . . is given as a series 
in the coupling, and each expansion coefficient is calculable in perturbation theory. Note that 
due to the fact that the same light-ray operators enter the definition of DA's and non-forward 



2 Or more general correlation functions of the elementary field operators in cases when the operator product 
expansion does not exist. 
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distributions, the corresponding evolution kernels are in one-to-one correspondence [T^. In leading 
order (LO) the calculations are quite straightforward since there exist different reliable formalisms 
for the evaluation of non-forward evolution kernels in the light-cone position as well as light-cone 
fraction representations. On the other hand the direct calculation beyond leading order turns out 
to be far from being a simple task which has only been solved for the evolution kernel in the non- 
singlet sector [19|, [2]J and was checked numerically in Ref . [22] . The most promising approach 
for two-loop calculations developed in Ref. J20j can be used for future explicit NLO computation 
in the singlet channel. However, the analytical solutions of the evolution equations for the DA 
IT , |24| , |25| , [26[| and non-forward distributions |16|, ^7\, [2S| are known in terms of the anomalous 



dimensions 7^ of the, so-called, twist-2 conformal composite operators. They can be obtained 
from the evolution kernels by forming moments with respect to Gegenbauer polynomials 

1 dxC u j {2x - l)V(x,y\a s ) -- 



~E7i*(«.)^(2y-i), 

Z k=0 



where the numerical value of the indices \i and v depends on the channel under consideration. 
Unfortunately, beyond LO this calculation is very cumbersome and could not be performed ana- 
lytically (only the first few entries of this matrix have been evaluated numerically |25], |23|). 

Owing to the difficulties we have sketched, every field-theoretical idea is welcome which can 
help us to face the above mentioned problems in a more economical way. On the first glance one 
could think that the field-theoretical background of QCD has already been well studied in the 
early days. However, after more than two decades of computation experience it is desirable to 
have a new look on certain issues. In this paper we apply an old idea of conformal symmetry^, 
which arises from the generalization of the Poincare group up to the fifteen parameter conformal 
group isomorphic to 50(4,2) — the largest possible group of transformations which leave light 
cone invariant. Besides translations and four-dimensional rotations it contains dilatations and 
special conformal transformations [^, |30|j . 

In the massless case the classical QCD action is invariant under infinitesimal conformal trans- 
formations. The restrictions coming from conformal symmetry in massless scalar and gauge field 
theories were studied already more than three decades ago when it was believed that this symme- 
try is broken only softly [f29. 30. 31. B3]. To give an example for the power of conformal symmetry, 
we remind the reader that in a conformal invariant theory the functional form of both elementary 
Green functions and current correlators is severely constrained: two- and three-point functions 
are fixed completely up to their normalization while n-point (n > 4) correlators are expressed in 



terms of a single function of anharmonic ratios [33 



3 Here we warn the reader from simple minded understanding of this terminology. From the consequent discussion 
it will be clear that our point of view is more sophisticated. It provides us with a powerful machinery, which enables 
us to predict certain quantities in a very simple manner. 
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Unfortunately, after quantization the theory is not conformal invariant even in the massless 



limit. The symmetry is broken due to anomalies related to the ultraviolet (UV) divergencies |34], ^5 



3(| |37], |38J which arise when one goes beyond the tree level. The regularization of these divergencies 
introduces a dimensional parameter which breaks conformal symmetry and this breaking still 
remains after renormalization and removing of the regularization. It shows up in the trace anomaly 



of the energy- momentum tensor |39], £4CJ, [41], [42| . Getting rid of the unphysical contributions in the 
trace (see below Eq. fl20|) ) such as BRST-exact operators, which originate from the gauge fixing 
and ghost Lagrangians, and equations of motion we end up finally with a symmetry breaking 
term proportional to the /3-function. If one assumes the existence of a non-perturbative fixed 
point g* with f3(g*) = 0, the theory would be conformal invariant. Formally we speak about 
the hypothetical conformal limit of the theory and set (3 equal to zero. Thus, the trace anomaly 
controls the conformal symmetry breaking in the Ward identities of the elementary fields operators 
|[35|| . However, this is by no means true when one is interested in more complicated Green functions, 
for instance, with composite operator insertion, since this one-to-one correspondence between both 
anomalies is lost. This means that in the conformal Ward identities (CWIs) the breaking of the 
special conformal symmetry is not controlled by the dilatation anomaly alone. Instead, new 
anomalies enter the game. 

This should sound a warning from the naive use of the conformal symmetry arguments. Indeed, 
in the past its application has led to a conflict between its predictions and explicit calculations 
beyond leading order ||43|| . For instance, in the above mentioned example of the renormalization 
of the conformal Ward identities with a composite operator insertion there arises a new symmetry 
breaking term induced by the UV divergencies occurring in the operator product of the composite 
operator and the trace anomaly. Such a term is not proportional to the /^-function. 

To study this breaking, we will construct Ward identities for the Green functions of the ele- 
mentary fields with conformal operator insertion which will serve as a main tool for the derivation 
of constraints for the anomalous dimensions beyond leading order. In general the renormaliza- 
tion matrix of the twist-2 conformal operators is triangular due to Lorentz invariance. In leading 
order it is diagonal so that conformally covariant operators do not mix under renormalization, 
i.e. 7^ = bjklf^ IQJ |24|, |44], Its eigenvalues 7j coincide with the Mellin moments of the 

DGLAP kernels. Therefore, the diagonal part[] of the exclusive evolution kernel can be obtained 
by the transformation [47| 

V D (x,y)= [ 1 dzP(z)f:^lq(2x-l)^C^2y-l), (3) 



4 The formula given below is very schematic and valid for direct (QQ or GG) channels. For the exhaustive 
treatment of this problem the reader is referred to the paper §47[ . 
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where Nj(u) = 2~ Au+l jS^^^t an d w {v\ v ) — {yy) v ~ 1 ^ 2 are the normalization and weight factors, 
respectively. Thus, conformal covariance at tree level is sufficient for the diagonality of the one- 
loop anomalous dimension matrix. However, beyond tree level conformal symmetry is broken 
even for (3 = 0. Of course, the scaling dimensions of the operators are changed by the anomalous 
ones but this shift does not ensure the covariance of the renormalized operators under special 
conformal transformation in a given renormalization scheme [|48"| , |49| , |50| , |47| . This fact is easy to 
understand, while the breaking of dilatation is determined in LO only by the leading logarithmic 
term, the breaking of special conformal symmetry is caused by both the logarithmic term and 
the remaining constant one. Since the logarithmic independent term is fixed by the normalization 
condition, which is implicitly given in the MS scheme, we can require a normalization which 
ensures conformal covariance of the renormalized operators^ |50| , [5TJ| . This procedure is equivalent 
to a finite renormalization and provides the, so-called, conformally covariant scheme. Therefore, 
the renormalized conformally covariant operators at LO possess a diagonal anomalous dimension 
matrix in two-loop approximation, up to a term proportional to the /3-function. Obviously, the 
non-diagonal anomalous dimension matrix in the MS scheme in NLO can be obtained by 0(a s ) 
finite renormalization of the two-loop forward anomalous dimensions. When translated into the 
language of evolution kernels the whole result reads now: 

V(x,y\a.) = V D (x,y\a s ) + V ND (x,y\a s ), V ND (x, y\a s ) = (g) V ND «(x, y) + 0(a»). (4) 

In other words the off-diagonal part \Z ND W of the kernel is a quasi-one-loop object and can be 
calculated in this order. This entry will be evaluated in the present paper. 

While in previous studies we restricted ourselves to the Abelian case, we extend in this paper 
the conformal machinery to the flavour singlet sector in QCD. Consequent presentation runs as 
follows. In the next section we summarize some basic formulae. In section 3 we derive the 



renormalized CWI analogous to the ones considered in Refs. [36 , p7| j3qj for the composite operators 



in the dimensionally regularized theory. We do not rely there on any assumption about symmetry 
breaking. The commutator constraints for special conformal and anomalous dimension matrices 
are established in section 4. Section 5 is devoted to calculation of the special conformal anomaly 
kernels in both light-cone position and fraction representations. Their transformation to the 
language of Gegenbauer moments is given in section 6. The subsequent section contains explicit 
formulae for two-loop anomalous dimensions as well as ER-BL kernels. Then we give the general 
formalism for the solution of the evolution equation in the required approximation. In the final 
section 8 we give our conclusions. To clarify the presentations given in the body of the paper 



5 Exploiting this fact together with the conformal operator product expansion we had, recently, the opportunity 
to give several predictions for some non- forward processes BQ, BU |32] . 
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we add three technical appendices. In the first one we list all Feynman rules required in the 
calculation of the conformal anomalies. In appendix B we sketch the light-cone position formalism 
for the computation of the leading order evolution kernels. We mainly focus our attention on the 
most complicated gluonic sector. The last appendix summarizes formulae for the hypergeometric 
functions which appear in the evaluation of the Gegenbauer moments of the LO evolution kernels. 



2 Preliminary remarks. 

To start with, let us recall some basic equations of QCD. The dimensionally regularized and 
renormalized Lagrangian is given by 

- f {f; v ) 2 - r hc F; v B\Bi - |w (f abc BlBtf - 1 faBi) 2 

+ ZsW + ^i^/By, (5) 

where d = 4 — 2e is the space-time dimension and /i is a mass parameter with dimension one 
introduced in order to keep the coupling constant dimensionless. We introduce here the symbol 
as a shorthand notation for the Abelian part of the full QCD field strength tensor G a ^ v = 
+ yfgX f abc B' \B„. The canonical dimension of the elementary fields are the following: d c ^ n = 
| — e for fermions, d C Q U = 1 — e for gluons, d^ n = d — 2 and d c ^ a = for the anti-ghost and 
ghost fields, respectively. This Lagrangian is invariant under the following renormalized BRST- 
t r ansf or mat ions : 

5 brst^ = -itf g z xU a t a i)b\ 5 BRST B a ^ = Z 3 D^u a 5X, 

^BRST^a = ^^gZJ^cu^X, 5 BRST u a = ^B^SX, (6) 

where SX is a renormalized Grassman variable. The covariant derivative is defined as follows 
Dfj, = — ifi e gXT a B^, where T a is the generator in the fundamental (T a 0j = or adjoint 

(T b <p a = if abc (j) c ) representations depending on the object it is acting on. The Ward-Takahashi 

_3 I 

identities imply the following relations between the renormalization constants Z\Z 3 2 = ZIZ Z = 
ZiZ 3 ~ 1 Z 3 2 = ZiZ 2 ~ 1 Z 3 2 . In the MS scheme Z-factors are defined as Laurent series in e: Z(g, e) = 

1 + Z^(g, e)/e + Z^(g, e)/e 2 + For our consequent discussion we take the renormalization 

constants Z 2 , Z 3 , Z 3 and X = Z\Z 3 X = Z\Z 3 X = Z\Z 2 X as independent ones. 

The renormalization group coefficients and the anomalous dimensions of the fields are 

Pe(g) = ^ = -eg + (3(g), a = A^lnf = -2j G , ^ = \^ lnZ ^ ( 7 ) 
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(<f) runs over parton species entering the Lagrangian {<fr = tp,G,u)}) where = Z2, Zq = 
Z 3 , Z Q = Z 3 . These coefficients are determined by the residues of the corresponding Z-factors 
]53[ and in lowest order approximation, | = f^(3 + Pi + . . ., 7^ = §^7^ + • • ., they are 
given by 

H <> =\T F N l - l lc A , 7 <°>(0 = \t„N, + f f ({ - H) , 7 f( f ) = |c F , (8) 

where the group theoretical factors are Ca = 3, (7p = 4/3, and Tp = 1/2 for SU C (3). 

In the following we are interested in the renormalization of the conformal operators in the 
flavour singlet sector, which we define as two-dimensional vector: 

°" - K) ■ (9) 

The group-theoretical construction of these local operators was given sometimes ago in Refs. 



44, [45| (see also Ref. |55] for a non-local version of conformal operators) 

= G +(1 (id + ) l - 1 {/""' )c]-A^r\o l - no) 



\Q V) 
QqA 



51 





where d=d + d and D=D — D- The + and — components are obtained by contraction with 
the two light-like vectors n and n*, such that n 2 = n* 2 = and nn* = 1. The indices of the 
Gegenbauer polynomials are determined by the canonical dimensions d™ n and the spins of the 
fields: v+ = df n + S<p -l/2. 

The composite operators (|T0| ) form an infinite representation of the collinear conformal algebra 
at the tree level. This subalgebra 577(1, 1) of the whole conformal algebra contains the generators 

V + = rvP M , M- + = n;M^n u , V, K- = n*/C M . (11) 

The total momentum V+ is a raising, the special conformal generator /C_ is a lowering and the 
dilatation T> and the Lorentz Ai-+ generators are diagonal operators: 

SlOn = i[O jh V + ]- = iO jl+1 , 6? + O jl = i[O jh M-+]- = -{l + l)O jh (12) 
6 s Oji = i[O jl ,V)- = -(l + 3)O jl , 5 c lO jl = i[O jh IC^=ta(jJ)0^ u (13) 

with ajk(l) = Sjk ■ a(j,l) and a(j, I) = 2(j — + 1 + 3). There is one Casimir operator of the 
collinear conformal algebra, which is given by J 2 = ^V+K,- — j(T> + M.-+) 2 — %(T> + AA-+). The 
operators introduced above possess conformal spin j + 1: J 2 Oji = (j + + 2)Oji, spin I + 1 
and scale dimension 1 + 3. 
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To study the conformal anomalies we consider the CWI for the connected Green functions^ 
defined by the formula 

nn m JD<j>[O jl ]<j>(x 1 )...<j>(x k )e iS 

{[ ° jl]X) = JD^S > ( 14 ) 

X is a product of elementary fields X = ni=i ^ipi) a ^ different space-time points. The integration 
is weighted with the factor e %s , where S is the QCD action determined in terms of the Lagrangian 
(|) S = J d d xC, and 

[O j i\ = Y,Z jk (e,g)0 1A (15) 

fc=0 

are the renormalized composite operators. Here we neglect possible gauge-variant counterterms, 
since they give no contribution to the physical sector of the theory. The four entries of the matrix 
AB Zjf.(e, g) with A,B = {Q,G} are triangular and independent of I due to Poincare invariance. 
In the MS scheme the expansion AB Zjf.(e, g) = 5 AB 5jk + AB Z^\(g) /e + ■ ■ • is valid [^, [54]. Since 
the regularization is manifest in the action, the true conformal Ward identities, i.e. those which 
are valid in the interacting theory, can be obtained by a formal conformal variation of Eq. ( |I4| ) : 

([O jt ]SX) = -(5[0 3l }X) - (i[O sl \5SX). (16) 

Here the LHS contains a differential operator acting on a renormalized Green function and, there- 
fore, it is finite by itself. Thus, the same property should be satisfied by the RHS. Because of 
Poincare invariance, the variations 5 B S and S vanish identically and thus, with Eq. (|i~2|) the 
corresponding Ward identities read: 

([O fl ]6*X) = -i([O jl+1 ]X) ([O^X) = {l + l)([O fl ]X). (17) 

In the case of conformal transformations, i.e. 5 = {5 s , 5^}, the variations of the action do not 
vanish anymore. The operator product [Ojj\5S appearing on the RHS of Eq. ( fI6|) is infinite. The 
removal of the UV divergencies induces the conformal anomalies, i.e. anomalous dimension and 
special conformal anomaly matrices. 

Before we start our study of the conformal variation of the action, we write down the variation 
of the composite operators when the interaction is switched on. Even in the interacting theory it 
is convenient to deal with scale dimensions which are equal to the canonical ones in 4-dimensional 
space-time, i.e. da — 1, — |. Therefore, the infinitesimal conformal transformations of the 



renormalized operators (]I5) are simply given by 

J \\Za(i)Z-\ iL 



5 s [O fl ] = -(/ + 3)[Cy, 5 C [O jt ] = i {^(O^ 1 },, Pki-i}- (11 



k=0 



6 It is straightforward to perform Legendre transformation to derive the CWI for the 1-particle irreducible Green 
functions. We will avoid this technical step, since we can obtain all desired results from the CWI of the connected 
Green functions. 
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The choice da — 1, d^ — |, d a — d — 2 and = ensures that the conformal variation of 
the action can be expressed by well classified operators: gauge invariant operators, BRST exact 
operators, and equation of motion (EOM) operators. A further advantage of this set of scale 
dimensions is that it allows to use Jackiw's conformal covariant transformation law for the gauge 
fields ||56|| . After some algebra, the final result can be written as 

5 s S = J d d x A(x), 5°S = J d d x 2x v A(x), (19) 

where the integrand A(x) can be expressed in terms of the regularized but unrenormalized trace 
anomaly of the energy- momentum tensor 0^ u : 

A(x) = -0^(x) - d G Q G (x) - d^Q^x) - d a Q Q {x) 

= t{0 A (x) + B (x) + Q Q (x)-Q^(x)} + (d-2)d^O B ^x). (20) 

Here we introduce the following set of operators 

7 2 JvBRST rBRST 

O a {x) = -± (G%) , O b (x) = —tfd^B*, Bll (x) = —UTS}, (21) 
as well as EOM operators 

It is desirable to express these variations of the action in terms of renormalized operators as 
this allows to neglect then terms proportional to e. For this reason the renormalization problem 



of the above mentioned operators have to be solved | 57| . The only gauge invariant operator 
Oa (class A) is of twist-4 and needs counterterms of class B operators, given as BRST-exact 
ones, and equation of motion operators (class C). The twist-4 class B operators are O b and the 
total derivative of the twist two operator d il P Bii . Since there is no other class B operator which 
they can mix with, only class C operators are needed as counterterms for the class B operators. 
This set of twist-4 and twist-2 operators is closed under renormalization. First, we discuss the 
renormalization of the twist-2 operator BfJi . Since there is no further class B or class C operator 
with dimension three and the same quantum numbers we may conclude that this operator is 
renormalized O b (x) = [O b (x)\. It has vanishing physical matrix elements at non-zero momentum 



transfer [58] while at zero momentum transfer it vanishes identically as it contains no poles. The 



renormalization problem of the twist-4 operators at zero momentum can be solved by a direct 



S 



calculation of the following differential vertex operator insertions^} 

[A*] = = l -{p B ] + a,,}, [a«] = g|-s = + \o B ] + n G + n B . (23) 

We tentatively expressed in these equations the derivatives of the renormalization constants in 
terms of the derivatives of independent set. From the finiteness of the differential vertex operator 

insertions the counterterms of [0^(2?)] and [(9b(x)] have been fixed up to total derivatives (here 

_i 

Z g = XZ 3 2 ): 

[Oa{x)\ = (l+g^^)^(x)+(g^-2e^){^(x) + ^(x) + f2 G (x)} (24) 
Pb{x)\ = (l + 2e^^)0B(x) + 2e^^{^x) + ^(x)} + e^^^(a:) (25) 



d\nZ 3 

Making use of the representation 



+2£^r^ Q Q (x) + dJl B (x). 



Ms) L 



d\nZ g \ 

+ g— ^-M (26) 



dg 

we can express e {O a (x) + O b (x)} in terms of the renormalized operators only 

e{0 A (x) + O b {x)} = -^[O a {x)\ ~ (j- 7g) [Ob(x)\ + (la - ^ {(2 G {x) + Q Q {x)} 

+ 7^, f2^(x) + 2juf2a,(x). (27) 

Thus, inserting our findings into Eq. ( p0| ) we obtain the renormalized anomaly: 

A(x) = -^p A {x)}~l^~ lG ^{[o B {x)]^n^x)}+l^ G -^n G {x) 

+( 7 ^ - e)fi^(z) + 2 7w /2 (z) + (d - 2)9 M [0 B/1 (z)]. (28) 



3 Renormalization of the Ward identities. 

Combining the results we have obtained so far in Eqs. ([IB]) , ([Of), Q19D and (p8|), we will come to 

the desired conformal Ward identities. However, since the products of composite operators possess 

7 Throughout the paper, to simplify notations, we use the following conventions: for unintegrated operator 
insertions and EOM operators we keep the dependence on the space-time point explicit <D(x); for the integrated 
quantities with weight function 1 we just omit this dependence O = J d d xO(x), while for weight 2x- we use 
Or =j d d x2x-0(x). 
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UV divergencies, we have to take care of their subtractive renormalization. The main problem 
here is that no theorems for the renormalization of (gauge invariant) operator products exist [j59fl . 

This section is devoted to the solution of this problem for the operator products appearing 
in the CWI (16). It is worth to mention that in our previous work on Abelian gauge theory ||7 



we discriminate a rigorous treatment in favour of a symmetrical handling of the quark and gluon 
sectors. Due to the non- Abelian character of the present consideration it is helpful to deal with 
the well classified operators Oi. Of course, both approaches — the one which we have used there 
and the other we describe below — lead to the same results for a theory with U{1) symmetry 
group. 

Fortunately, due to the fact that the products involved are composed of gauge invariant con- 
formal operators and differential vertex operator insertions the form of the gauge- variant part can 



be fixed from the study of the Green function (14) with the help of the action principle. So that 
finally we have for the unintegrated operators: 

i[0 A {x)][O fl ] = i[0 A {x)O fl ] - £ {z A }. k [O m ] - £ {Zl} jk Pu-i] 



k=0 k=0 



( <91nX d\nX\ nl , 5 T/n n 



;[o B .M][cy = iiOrWOj,] - s<-0{x) Y. {z B } , [o u ] - '-a + s^(x) £ {zi}.. [o*,-,] - . . . 

fc=0 3 / fc=0 

-(4")^^ (30) 

The ellipses stand for the higher derivatives counterterms which are not of relevance for our present 
study. For convenience we introduced here the operator combination [(9b. (x)] = [Ob (a;)] +Q Q {x), 
which on tree level is the gauge fixing term and thus require only counterterms proportional to 
£. In the Landau gauge the operator product [Ob*(o;)][0^] is finite, since all counterterms will 
vanish. Obviously, the term containing the variation of the composite operators with respect to the 
gauge field is gauge- variant (it cannot be decomposed in terms of A, B or EOM operators). In the 
following subsection we show that the appearance of this term ensures the correct renormalization 
group equation for the composite operators. We also will fix the Z-matrices and Zg in terms 
of known renormalization group coefficients. 

Fortunately, by employing BRST invariance and naive power counting we are able to show in 
general that the operator product [Oji] [Ob^x)] is renormalized: 

[O jl ][0 Blt (x)) = [O jl O Bll (x)]. (31) 

Finally, we have also to deal with the product of composite operators and EOM operators. 
For these products we do not use the MS prescription, rather we perform a mere integration in 
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the path integral which results in equation: 

i[O sl \n^x) = i[0^(x)] - ^x)-^—^}. (32) 

where it is obvious that 

([O fl fy(x)]X) = iftOMW-^X). (33) 

For instance, for the counterterms of the products containing the EOM operator of the quark 
fields we have 



j 



(34) 



/ ddx2x - V ~k + = 2i t{^PQZ- l }. k [O kl -,\. (35) 



fc=0 



We have used above the projectors on the quark and gluon sectors 



,0 0/' V0 1 

and the 6-matrix with elements given by 



, 1 0\ P /0 0\ 

P Q = L n > P g = L n I (36) 



bjkM = 2 I" dx^l^ k (2x - 1) 
Jo Nu(b>) 



C v A2x - 1) 



/o ' N k (u) 

= 2 9 jk {(l + k + 2u)5 jk - [I + (-l) j ~ k ](k + v)} . (37) 

Due to the definition of the conformal operators (|T0|), the same matrix bjk(l) = bjk(l, 3/2) = 
bj-ik-i(l — 1)5/2) appears in the quark as well as in the gluon channels. 

3.1 Scale Ward identity. 

The renormalization of the dilatation Ward identity is quite straightforward in terms of the dif- 
ferential vertex operator insertions ([23|). Nevertheless, we present here an alternative way, which, 
however, can also be used for the renormalization of the special conformal Ward identities. After- 
wards we show that the appearing dilatation anomaly is just the anomalous dimension matrix of 
the composite operators. 

Taking into account the subtractive renormalization procedures in Eqs. (p9|), (|30|) and fl32|) , 



we find with Eqs. (|20| ) and (|2~8|) the following equation for the operator product appearing in the 
Ward identity (|T6|) (with 5 = 5 s ) expressed via renormalized insertions: 



i[0 3l }5 s S = j2{-ZA+(--l G ]z B -2( llp -e)ZP Q Z- 1 \ [O kl ] 
k=o I S \ S / J jk 

-^[OjtOA] - (| - 7G ) i[O fl OB.] + ( 7G - fj i[O fl n G ] (38) 
+( 7 ^ - e)ip 3l Q^} + 2 1 J[O jl n^\. 
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Inserting this result in the Ward identity fll6| ) we get with Eq. ([13]) in the limit e — > 0: 



{[O^X) = E{(' + 3 )l+7}., ([O kl ]X) + ti(i[O jl O A )X)+ (^-jg \{i[O jl O B ,]X) 

k=0 3 & \S J 



-^(iiOj^X) - (j G - (ip^X) - 2 lw (z[O jl Q ai }X), 
where the finiteness of the LHS of the Ward identity ensures the existence of the matrix 



7 



hm {-jZa - (f - TgJ Z b + 2( 7 ^ - ^ZPqZ' 1 ^ . 



(39) 



(40) 



It remains to show that the dilatation Ward identity (|39|) is nothing else than the familiar 
Callan-Symanzik equation [^0], |6l| and that 7 is the anomalous dimension matrix defined in the 
MS scheme. By differentiation of the renormalized Green function with respect to the parameters 
g and £ we define the differential vertex operator insertions [C^A*], so that 



3 ( r)7 

i[O iI ][A«]=i[O iJ A*] - Y\£-qZ Z ~ X 



k=0 
+ (l+g 



dZ_ 
dg' 

dhiX" 



dg 



3 ( f)7 

i[o Jl m= i [o Jl AZ] eU-^ 1 

k=o 1 a< s 



+ H 8 -^fz f 1 / d*xB?.(x) 



5B"(x) 



2P G Z 



- 1 } Pm]. (42) 

J jk 



We learn from Eqs. ( 23], 29, 30,41,42]) that the renormalization constants Z; t can be expressed in 
terms of the renormalization matrix of conformal operators and read 



dg 



dg 



zP G z-\ 



zP G z-\ 



(43) 
(44) 



Inserting these findings into the renormalized dilatation Ward identity (|39"D provides us the 
well-known form of this Ward identity: 



{[Oji}5 s X) = £{(/ + 3)l + 7 }. APki]X) + ^(ip^jX) + aiip^X) 



k=0 



and from Eq. 



-^(iio^n^x) - lG {i[o fl nG\x) - 2 lu} (i[o jl n a \x), 

follows the familiar definition for the anomalous dimension matrix: 



7 = 7 Z + 27 G P G + 2^P Q , with ^ z = -fi—ZZ' 1 = g—Z^. 



(45) 



(46) 
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3.2 Special conformal Ward identity. 

The renormalization of the special conformal Ward identity can now be performed along the 
same line as in the previous subsection for the dilatation. For the renormalization of the product 
[Oji)5^S we obtain from Eqs. (|B5p, fl3§ together with fl^, (§g), §H§ and (0): 



k=0 



g 



g 



^-jG)i[O jl O^] + ( lG -^]i[O jl na] 



g \ 6 / \ 6 / 

+ ( 7 ^ - e)i[0 3l n~^} + 2 7w i[O i ,/2-] + (d - 2)i[^A- t ], 



(47) 



where [A^ t ] = / d d x 2x_<9 / JO£ M (x)]. Inserting this equation together with the special conformal 
variation of the composite operator, given in Eq. (|18D , into the Ward identity (fl6|) (with 5 = S_) 
we obtain the desired CWI 



J 



'P 



k=0 g \g ) 



2(i[O jl A;jX) - ^(ip^X) -\ig--J (iPn^X) - 2 1 MO fl n-\X). (48 



The anomaly appearing here is called special conformal anomaly matrix and is defined by 



lim < Z 



a{l)-2( 7i) -e)P Q b 



Z- 1 - ^-Z A -(^- lG ) Z B - a(l) 



-2 llp P Q b + 2 



z^KPob 



1 - 



Mll- 



(49) 



where the second equality holds due to the fmiteness of 7 C . Note that this definition looks similar 
to the anomalous dimension matrix defined in Eq. 



4 Commutator constraints. 



A priori there is no way to relate the new renormalization matrices Z A and Z% to the renormal- 
ization group constants. Obviously, the special conformal anomaly matrix 7 C (/) contains therefore 
new information. However, from the conformal algebra it is possible to derive the following con- 
straints for the conformal anomaly matrices: 



[JC-,V+]- = -2i(V + M- + 
[£>,£_]_ =iK. 



f{l + 1) - 7 C (0 = -2 (7 
7,a(0+7 c (0 + 2-S(/)P ( 



2^; 
g j 



0. 



(50) 
(51) 
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The constraint ( pip tells us that the breaking of translation invariance by the infinitesimal special 
conformal transformation is controlled by the commutator relation between /C_ and V+ providing 



the spin dependence of the special conformal anomaly matrix. The constraint (51) relates the 
off-diagonal matrix elements of 7 to the matrix 7 C (Z) + 2^b(l)pQ. Since a(l) is a diagonal matrix, 
which is independent of the coupling a S} the n-loop approximation of 7 ND is determined by the 
{n — l)-loop approximation of 7 C (/) and /3-function. 



4.1 Constraint for the special conformal anomaly matrix. 



From the definitions (p~2|JT~3|) it follows that the commutator [/C_,"P + ]_ = —2i(T> + -M_ + ) corre- 
sponds to the following relation for the variations = 2(5 S + S M + ). We will employ this 
identity in order to prove the constraint (^0) which comes from the conformal algebra and the 
Ward identities (0) and ©: 



([0,,] [5^5 c ]_X) = ([O jl ]2(5 s + 5 M + ) X) 



£ {<l + 2)1 + 2 7 } ([O kl ]X) + 2^(i[O jl O A )X) + 2 ( ^ - lG } (i[0 3l O B «}X) 



P 



P 



^{ipflQ^x) - 2 [j G - ^ {i[o jt n G ]x) - 2 lu (i[o jl n ]x), 



(52) 



With the help of the Ward identities (|17]) and (fS|) we find immediately for the first term in the 
commutator on the LHS of Eq. (|52|): 

([O 0l ]5 p + 5 c X) = - X {a(/ + 1) + 7 C (/ + l)} jh ([O k i)X) + §([O iI+ iOl]AT> 



fc=0 



'P 



+ (j ~ 7gJ ([O il+ iOi.]^> - 2([0 j7 A cxt ]*> - 7^+1%]*) 



7g 



([O i H. 1 ^]^)-27 w <[O i ,+i^]Ar). 



(53) 



While for the second term we get 



P 



([0 3l ]6 c 6^X) = - X {&(/) +7 c (0) jfc ([0«+i]AT> + £<i[0*0A]OO 

fe=0 o 



lG - £\ iip^SlX) - 27 U <«[CW]<?*>- 



(54) 



To proceed further we have to get rid of the variation sign on the field monomial, i.e. 5+X. 
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Using the definition of the subtractive renormalization fl29|), fl30|) for the product of two com- 
posite operators we have 



-{iSlPflO^X) for i = {A,B*} 
([O jl+1 0-]X) - 2{i[O jl O i ]X) - ^ lA P G l\0 3l \X). 



(55) 



Note that the last term arises due to the non-minimal nature of the subtraction: the renormal- 
ization matrix Za consists of a finite term apart from divergent contributions oc -. 
Since the operator [A^ t ] = — 2n* J d d x [Ob^ is translation invariant we can write 



m^-js^x) = ([o jl+1 A:jx). 

For the terms containing the EOM we find immediately 

(i[Ojini]5^X) = 2([Oji] I dx(f)(x)-r-^-—X) - {ip^Sl [ dx2x^(x) T ^ nx X) 



x 



S(j)(x) 



([o jl+1 n;]x)-2{i[o jl n 4> ]x). 



(56) 



(57) 



Inserting our findings (|53f 56) in equation (52) and using the definition of the elements of the 
a-matrix yields 



E \ni + 1) - m + 2 [i - 2-Pg) ) ([o kl+ i\x) = o. 

k=o I V 8 I ) jk 



(5* 



Owing to linear independence of the Gegenbauer polynomials we can safely omit the Green func- 
tion and get finally the matrix equation (|50|). 



4.2 Constraint for the off-diagonal anomalous dimension matrix. 

Now we prove the commutator constraint (|5lD. Since 5 s = —fi-§-, instead of dealing with the 
corresponding commutator we apply the renormalization group operator 

d « 

D = fi— 1 + 7 
djj, 



on both sides of the special CWI (|48| ) without taking into account the differentiation of the field 
monomial X with respect to fi since these terms vanish by means of the special conformal Ward 
identity. The LHS vanishes identically since D[Ojk] = 0. To calculate the RHS we employ formally 
the renormalization group equations of the appearing operator insertions. For this purpose we 
need the following commutators 



d d 
^ dji ' dg 



d p e \ d ( d \ a 

: x 8^~ + 2 g^-lc fe, 

dg g J dg \ dg d£ 



dfi ' <9£ 



2 Nh) 4 (59) 
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Using them properly when one acts with the differential operator //j- on the operator insertions 
[O^], we obtain the desired renormalization group equations: 



g 



d_jT 
dgg, 



d , d 
2£ 



dg 



-la )[Ob*] 



(60) 



(61) 



where 







Now we list the action of the operator D on every term appearing on the RHS of the special 
conformal Ward identity. The first equation is obvious: 



k=0 



k=0 

d 



j, a + f + 2-P Q h 
g 



(62) 



+ Pi- 4 1] " + Z l B l ~ - 2^P Q b \ ([O kl ^)X). 



dg 



Here we have taken into account that the special conformal anomaly matrix is a gauge independent 
quantity. Next we calculate the action of D on the renormalized operator products by using 
the subtractive renormalization prescriptions (|29|) and (^) as well as the renormalization group 
equations derived above: 



D^([O fl 02]X) 
(3(8 



(63) 



g h - 2 4) ( 7g - -g) {[0ji0 ~ B * ]x) + f ? (4 ~ 2 4) w 



p4 



o 



[i]- 



<9g d£ 



D 



'P 



- la) (P !t O- B .}X) (64) 



+ 



'/3 



-7G E 

/ fc=0 



■d-g ZB , 



Due to Eq. fl33|) the following relation holds 



d 



D ^(g)([0 3l ^]X) = U-^( g ) ([<W]*>. 



(65) 
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In the above equations we have omitted the divergent contributions on the RHS (which cancel 
each other) since the LHS is finite. All other terms in the special conformal Ward identity vanish 
when the .D-operator acts on them. 

Collecting all terms listed above, we find the desired constraint 



k=0 



7,a(/)+7 c (0+2^(/)P Q 



> ([O kl - 1 ]X) = 0, (66) 

jk 



which implies finally the matrix equation (plf). 

We have also found the following expression for the gauge-dependent renormalization constant 

= 4 (4 7 *) ^ (6?) 

5 Leading order conformal anomalies. 

In this section we compute the conformal anomalies to LO in a general covariant gauge. Since we 
are interested in the present paper in the NLO approximation for the anomalous dimension matrix 
of conformal operators, we restrict ourselves to the LO approximation of the special conformal 
anomaly which reads 



r 



2 [wz w , h] _ + QQ z [ l ] ~ -2b QG Z t 



Here we have taken into account that in one-loop approximation the following equalities hold 

QQ Z$- = 2 lf b and QG Z l } ] - = 0. (69) 

For technical reason we perform the calculation in the light-cone position formalism. Then 
it is quite straightforward to transform these results to any other representation. We define the 
relevant light-ray operators as follows: 

0(ki,k 2 ) = , (70) 



where 



\Q Q v] 
\Q A 



\ G o 

Gq 



V 



g O(k u k 2 ) 



[ki,k 2 ) = ip(K 2 n)l + > §{n 2 n, Kin}ip(Kin) =p ip(nin)\ ^ + \ $[/«in, K 2 n]ijj(K 2 ri), (71) 
1 7+75 J 17+75] 

, k Ki,k 2 ) = G +ll (n 2 n){ ^ ) $[« 2 n, Kin)G v+ (Kiri). (72) 



le Liv — hj 
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A 



A 



(a) 



(b) 




Figure 1: One-loop graphs contributing to the renormalization constant Z A . 



The path-ordered link factor $[x2,xi] ensures gauge invariance and reads 



1 



drT a Bl{x{r))^ 




JO 



(73) 



where the path is parametrized as x(t) = rx 2 + rx\. The local conformal composite operators 
(|ToD are obtained by differentiation with respect to K\ and K2, namely: 



where d + = d Kl + d K2 and d + = d Kl - d K2 . 

In the next subsection we derive the special conformal anomaly matrix for the quark sector 
including chiral-odd (transversity) twist-2 operators. This matrix turns out to be universal in 
a sense that it is the same for all Dirac structures. Then we comment shortly our results for 
the mixed QG and GQ channels, which were derived previously in Abelian theory. In the last 
subsection we evaluate the special conformal anomaly matrix for the GG channel for parity even 
and odd operators. Again in all channels we observe the universality of these kernels, which implies 
that the additional special conformal symmetry breaking is the same in vector and axial-vector 
channels: the difference arises from the different scale anomalies, i.e. the one-loop anomalous 
dimensions of composite operators. 

5.1 Once more about quarks: chiral-even and -odd sectors. 

Although the quark sector was considered by us in detail in previous studies, nevertheless, it is 
worthwhile to address this issue from our present point of view. Moreover, we complete the list of 
two-loop anomalous dimensions by adding those of the chiral-odd (transversity) sector. Thus we 
should evaluate the renormalization constant from diagrams in Fig. |l|. This calculation is 




(74) 
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quite trivial. Using the Feynman rules given in appendix |A] we get easi 



a, % 



dz I dy 
(a) 2n e Jo Jo 



l[0 A ][ Q O r {K U K 2 )} 

x2 QQ )C( a) (y, z\S = Q) (d k J + d k j) [ Q O r (^i + yK2 , ZKl + zk 2 ) 



(75) 



l[0- A ][ Q O r {K U K 2 )] 



a. i 



(b+c) 2tt e Jo 



dz / dy (76) 
Jo 

2 QQ !C( b+c) (y, z\Z = Q) (<9 fcl + d k j) + QQ IC w (y, z^O^y^ + V k 2i zk x + zk 2 )\. 

The separate contributions Q®K.[ to the evolution kernel correspond to the Landau gauge, so that 
the sum Q®JC F = ®®K£ a \ + QQ /C^ +C ^ is gauge invariant as it should be, since for £ = the one-loop 
quark anomalous dimension vanishes, 7^°' ) (£ = 0) = 0. The /C^-kernel reads 



QQ )C w (y, 



-2C f 



ki- 



8{z) + — 
k 2+ 



(77) 



where we have introduced the plus-prescription for a singularities l/y k (with k = 1, 2, . . .) 



1 



fc-i 



Z 



k-j ' 



(7S 



When Eq. (|77| ) is transformed into the light-cone fraction representation (see Ref. [p2jl ) we repro- 
duce the known expression for the w{x, y)-function 



3 w(x, y) 



-2CV 



x — x) fx — > 
y(x-y) 2 \y^y / 



(79) 



This kernel is the same for all quark operators irrespective of their Dirac structure: J 1 = 7 + , 7+75 
or o+j_75. Therefore, the special conformal anomaly matrix reads < 3 < 3y c (°) = —faQOyi ) 4- QQy). 



5.2 Mixed channels. 

Note, that in the present approach there is no need to perform any calculation for the mixing QG 
channel. The u>-kernel for the GQ-channel was calculated by us in Ref. |I7| with the same result 
for the parity even and odd channels 



a *K? (y,z) = C F [5(z)-5{y)]. 
The momentum space kernel reads 



(80) 



1 



x — ► x 



GQ w(x,y)=2C F (-e(y-x) , 

\y \y-+y. 



Note the extra factor of 2 in front of the light-cone fraction kernel which comes as a result of the 
transformation from the coordinate to the momentum space according to Eqs. QA.l , A..2Q . 



3 The last argument of the kernel stands for the value of the gauge fixing parameter. 
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(/) (g) (h) (0 V) 




(k) (I) (m) 

Figure 2: One-loop diagrams which give rise to divergencies in the product of the renormalized 
operator insertions zfO^jpC^]. The symmetry factors s a for the diagrams are the following: 
— 1 for ex — a, b, c, /, g, h, i and s a = \ for a = d, e,j, k, I. The diagram (m) is identically zero 
in dimensional regularization. 

5.3 Gluonic sector. 

Using the representation of the conformal operators in terms of the derivatives of non-local string 
operators we can write the divergent part of the operator product i[(9^ 4 ][ G '(9(fi;i, ^2)] via the fol- 
lowing relation in leading order of the coupling constant: 

ip-^OiK^Ki)} = -L- dz dyl^Aiy^TOiyKt + y^zKx + zKz)} 

2tt 6 Jo Jo I 

+ GG iCA(y,z)Jd d x2x^(x)jJ^[ G 0(yK 1 +yK 2 ,z (82) 

Here the kernels IC(y, z) are the subjects of the calculations. The diagrammatic representation of 
the corresponding graphs is given in Fig. |2| and the corresponding Feynman rules can be found in 
Appendix |A]. The details of the calculations together with the general approach to the evolution 
equations for non-local string operators are given in Appendix 0. 
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Note that the factor in front of the operator Bj^[0] is gauge independent in leading order. It 
can be conveniently represented as 

1 a s [fa x 



dlnX d\nX 
-26 



dg 



e27r 



The result of the calculations given in detail in Appendix M can be summarized as follows 



-z2(K! + n 2 ) GG IC(y, z) - z(«! + « 2 )/3 %)^) + GG IC w (y, z) 
i(y-7g } (e = 0)) 



where the parity-even and -odd gluonic kernels read |63], |55], [16, ^4 

1 



C r ^4(l- 2 /-2 + 3^) + | 2 /-2 + 

+ <z-2 + 



S(z) 



0o 



GGtrA 



K A (y,z) 



GGtrV 



'JC v (y,z)-12C A yz, 



3) 

34) 
35) 

36) 
37) 



and 



GGirw 



% w (y,z) 



Ca 





V 




■ 1 " 




T 




" 1 " 








- 2 






- 2 




.) 




_z_ 


+ 




.2/. 


+ 


y 2 . 





*(*), 



which is the same in both channels. The light-cone fraction function reads from this 

x 2 9(y — x) 



GG w(x,y) = -2C A 
Using the identity 



y 2 (x - y)' 



1 -9(y-x) + -S(y-x) + 



V 



v 



y^y, 



(89) 



{id + ) l C»( d+/d + ) {{ Kl + K2 )0( Kl , K 2 )} 



Ki=0 



(90) 



k=0 



with the fe-matrix being given by Eq. (|37|), we can represent the above result (|4]) in the usual 
matrix form 

GO,, (91) 



GGf^c(0) = _gC(^(0) +A) g + „ Mu , i 



where the w-matrix possesses elements defined as Gegenbauer moments of the kernel (BW). This 
quantity will be evaluated in the next section. Note also that Eq. ( pij) as well as other conformal 
anomalies fulfill the constraint equality (|50|) . 



6 Gegenbauer moments of the evolution kernels. 

This section is mainly technical so that rather experienced reader can skip this presentation. Here 
we intend to give the method for evaluation of the moments of the evolution kernels we have 
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derived in the preceding section in the basis of Gegenbauer polynomials. We deal here only with 
gluonic kernels. 

Let us demonstrate the main technical steps for the integration of the functions with +- 
prescription, e.g. the first term in Eq. (|39"D, GG w^(x, y) = [x 2 /y 2 9(y—x)/(x—y) 2 ] + +(x^x,y^y) 
. Integrating by parts over x we decrease the power of the singularity and end up with equation 
(according to Eq. (|C.3|) we have replaced the Gegenbauer polynomial by 2 P*i) 

T(j + A)T(k + A) /•!_,__ w(y\^ 



1+ -1 



1 r 2 (5)r(j)r(A;) J *iv fc _ U2 

where the integrals Jj read 
JM = 



-k+l,k + A 
3 



yj{J 1 (y) + 2J 2 (y)}, 
(92) 



dx 



x 



1 — X 



d 



d(xy) 



My) 



dxx 



X 


1 ■ 


.1 — X. 


+ xy . 



2 -Pi 
2F1 



xy 

xy ) - 2F1 



-3 + 1,3 + a 
3 







(93) 



-3 + 1,3+4 
3 

-3 + 1,3+4 
3 

Since at the end we get the result in terms of hypergeometric functions and their derivatives with 
respect to the lower index, the simplicity of reduction to elementary functions depends on the 
handling of the integrals Jj. Therefore, at this step it is instructive to expand the integrand of Jj 
with respect to a complete set of Gegenbauer polynomials. Namely, 



d F 

dx 



1 

X 



-3 + 1,3+4 
3 

-3 + 1,3+4 
3 



x 



X 



e n 2^1 



1=1 
2F1 



-1 + 1,1 + 4 
3 

3 + 1,3 + 4 
3 



X 







1=1 



The integrals (|93|) can be easily evaluated with the result 



OO Q 

1=1 ut 



2 -Pi 

e=0 



-I + 1,1 + A 
3 + e 



2 -Pi 

e=0 



-I + 1,1 + A 
3 



-I + 1,1 + A 
3 + e 



where we have used analytical regularization of the singular distribution 



Thus we get finally for GG w^ 1} 



[1-x 

jk 



dx F{x) = g 



dx 



(1 — x) 



— {T{x)-T{l)}. 



[1 + 



\j-k] 



{fji h ik(j) - ( 2 fji + Cji)gik(j)} 

1=1 



(94) 
(95) 

(96) 
(97) 

(98) 



with 



hik{j) 

gikU) 



r(j+4)r(A; + 4) /-i w(y\ 



r*(5)T(j)T(k) Jo "N^m 



dy- 



^2 



Fi 



-k + l,k + 4 
3 

f 1 



x 



d_ 



t=0 



-1 + 1,1 + 4 
3 + e 



y .(99) 
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Now these matrices can be easily computed. The main steps are reduced to the use of the repre- 
sentation ( |C.4j ) and integration by parts until all derivatives will act on the right hypergeometric 
function. Then with the help of Eqs. ( U.5 ) and (|C.2|) , we obtain the result in terms of hypergeo- 
metric functions of argument 1. The reduction to elementary functions is accomplished according 
to the formulae ( |CJ . (j| - |CJ . lOj ) . In this way we get 

r(j)r(/ + 4) 



fji 
hik(j) 

9ik(j) 



-[l-(_l)J-%_ lfl (3 + 20 



r(j + 4)r(/)' 



l u ( iv-ifl (o i on r(j)r(l + 4) 
-2 [1 ~ ( - 1} ]^-M(3 + 2/) r0 . + 4)r(/) 



r(j + 4)r(p 
r(j)r(/ + 4) 



s, 



Ik 



r(j + 4)r(p 
r(j)r(/ + 4)' 
r( J + 4)r(0 



(100) 
(101) 
(102) 
(103) 



Combining these results together we have 

GG w$ = [1 + 0^(3 + 2k){{^{2) - i/j(k + 2)) 

(3 + 21) 



- 1 



r(j + 4)r(fc) 
row + 4) 

T(j + A)T(k) 



+ y [i - (-i)^ fc , 

4 1 J J (/-A;)(/ + A; + 3) \T(j)l { k : 1) 
The final summation can be easily performed with help of the formulae 

(3 + 21) 



+ 1 



(104) 



J (/-^)(/ + fc + 3) 



«=/t+i 



•i 



(3 + 2/) 



= A^ + V0 + 2)-V(A: + 2) 

r(i) r(*) 



(105) 



(/- fc)(Z + fc + 3)r(/ + 4) T(A; + 4) 



A ife - VO + 2) + + 2) 



0-fc)Q + fc + 3) 

+ 1)0 + 2) 



(106) 



Adding the contributions from the last two terms in Eq. (p9 



<«>..^ - [1 + (-l)^]^._ 2ife (3 + 2k)\-^_{j - k)(j + k + 3) [0 - k)(j + k + 3) + 2](107) 



2T(A; + 4) 



GG w f) 



jk - -[l + (-iy-^._ 2jfc (3 + 2A:)- 



rp + 4)r(fc) 
row + 4) 



(108) 



we get the quantity in question 

-2C A [l + (~iy- k }9^ k (3 + 2k) 

~r{j + i)r{k) 1 



GG r 



2Ca£- 



(109) 



x { 2A jfc + (A jfc - VO + 2) + V(l)) 



row + 4) 



+ 20-A;)0' + A; + 3) 



V{k) 



T(k + A) 
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where we have introduced the matrix A the elements of which are defined by 



.4 



'j + k + 4 N 



'j - k" 



+ 2^(j-A:)-V0' + 2)-V(l). 



;no) 



The matrices for the QQ and GQ cases can be found analogously and were derived in our 
previous studies [P 



-2C f 



l + (-l)i-*]0 i _ 2)fc (3 + 2A:) 



x {^ + (^-^ + 2) + ^(l)) ^-^ ( t^ 3) 



-2CV 



l + f-l 



%_ 2ifc (3 + 2fc) 



l (j-fc)(j + fc + 3) 
6 (A; + l)(fc + 2) 



(111) 
(112) 



In addition ® G w = 0. 



7 Two-loop anomalous dimensions and evolution kernels. 

From the consistency relation (|50f) we can now derive the equation which defines the off-diagonal 
elements of the two-loop anomalous dimension matrix via the one-loop quantities evaluated in the 
preceding sections. Namely, we get (j > k) 

a {J^njk —lj Ijk ~^jk Ik +P0\1j ^Q-^Qlk )°jk- l-L-L^J 

The complete entry which governs the evolution in NLO is defined by the sum 7-^ = ^fjSjh + 
T^P , with 7^ being the two-loop forward anomalous dimensions 0, |8|, U • Defining new matrices 
djk = bjk/a(j, k) and = Wjk/a(j, k), we can rewrite the above equality in components 

= (^f - <^») (ft -^™) (114) 



, QG (0)GG _QQ QG (0) 

+ Ik 9jk 9jk Ik > 

GO ND(1) _ /GO (0) _ GO (0)\ , ( R QO,(0)\ _ f GgJV _ GGj°)\ j GO,(°) 



i-r = - UQ <) ^ (#> - w 7rj - [ GG i7 - uu ik') ik (He) 

+ GQ if )Q % k - G % k GQ if + (Sf - QQ ik 0) ) G % k , 

GG iT ] = faf-^Mfh-^+^gJ (117) 



GQ lf- GQ lf)d 3 k QG lf- G %k QG lf. 
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Here the leading order anomalous dimensions of conformal operators read [|65], [56], [J7|, i7| 

2 



qqJ°) 

i 3 



3 



1 3 



GO r (0) 



-Cp 3 + 



(i + l)(i + 2) 



-2AN f T F 



x 



-#0' + 2) + 4^(l)J 

/ J 2 + 3j + 4, for even parity 
j(j + l)(j + 2)(j + 3) " \ j(j + 3), for odd parity 

—Cf j j 2 + 3j + 4, for even parity 

3(j + l)(j + 2) X j j(J + 3), for odd parity 

M 

8C A 



-C A -4</>(j + 2) + 4^(1) 



(118) 
(119) 
(120) 
(121) 



j J + 3j + 3, for even parity 
"i0' + l)0' + 2)(j + 3) X j i(j + 3), for odd parity 



and for the transversity sector |5T|, |68 



-O(3-4^0' + 2) + 4V(l)). 



122) 



Note that Eqs. ( TOM) fulfill the Dokshitzer SUSY relation @, g (here Ca = Cf = 27» 

QO/) + = 3_QG (°) , GG (0) 



(123) 



The prefactors j/6 (6/j) result from the conventional normalization of the Gegenbauer polynomi- 
als. 

When translated into the language of the evolution kernels their off-diagonal parts in NLO 
read: 



OOyND(l) = _ (J _ v) \QQy ^ (QQyiO) + fhj\ 

+ QQ g g, _ g QQ^ + QGp, ^ GQ y (0) _ QGyW ^ GQ g \ 

-(I-V)\® G V® ( GG v (0) + ^i) 

+ QQy g QG^(O) + QQ^ QGy(O) _ QGy(O) ^ GG g 

-(i-v){ GQ v® (W (0) + ^/) 



QGy 



GQyND(l) 



(124) 



(125) 



GGyND(l) 



(126) 



+ GG^ GQ^(O) _ GQ^(O) ^ QQ^ + GG^ ^ GQ^(O) _ GGy(0) ^ GQg + GQ^ ^ QQ^(O) 



(X-D)( ( 



GG U (g) GG 



F (0 + |/ 



(127) 



+ GG g (X G< ' 



y(0) _ GGyiO) ^ GG + GQy ^ QGy(°) + GQ ^ QGy 



(0) 
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Here we use a shorthand ® = Jq 1 dz so that for any two test functions Ti®T2 = J 1 dzT\{x, z)t?.{z, y). 
As usual, the (X — D)-projector extracts the off-diagonal matrix elements of any function in the 
basis of the Gegenbauer polynomials ||7 . 

Let us mention that the 75-problem which might arise in dimensional regularization does not 
show up in the off-diagonal elements in NLO. This follows from the fact that the finite renormal- 
ization constant which restore the chiral invariance broken in the t' Hooft-Veltman-Breitenlohner- 
Maison scheme is diagonal in the basis of the Gegenbauer polynomials and thus cannot affect the 
results given above |j9 . 



The normalization of the Gegenbauer moments is defined as follows 

f 1 dxCf A \2x-l) AB K(x,y) = j: AB K ]k Cl {B \2y-l\ 

JO t^n 



;i28) 



k=0 



so that we have the correspondence K = V and Kjk = —oljk, K = (J — V)g and Kj}. = Qj-2,kQjk 
K = (l-V)V and K jk 
following form 



@j~2,k( AB jf* ) — AB lk^)^jk- The leading order ER-BL kernels have the 



ab t ^(o). 



X — > X 



l V^(x,y)=9(y-x) AB F(x,y)± 

withP] the following kernels for the parity-even [HJ [LI , |65|, 1(| [15 
B3, K7H channels 



for 



A = B 



(129) 

47j and -odd |7], [TfJ [L6|, p 



QQ F(x, y) 



x 



C F - 

y 



1 + T^u + \ i(x - y) 



nr^, \ , T m x (2x — y — l, for even parity 

QG F(x,y) = 2N f T F —\ _ P y 

1 —y, tor odd parity 



GQ F(x,y) 
GG F(x,y) 



C F - 



y*y 

x \ 2y — x, for even parity 



y I x, 



C A 



x 



for odd parity 

1 A) 



(130) 
(131) 
(132) 



(y - x)- 



2C, 



-5(x - y) 



+ 2C A 



x 2 ( x + y(l + 2x), for even parity 



JJ 



While for the chiral odd sector we have (see also 



<*>F{x,y)=C F - 

y 



7 — — x — o 5 i x ~y) 

{y-x) + 2 



for odd parity 

(133) 



(134) 



For the dotted kernels we get E7 
AB V(x,y) = 6{y-x) ^ w 



x 



Fix, y) In - + A AB F(x, y) 

y 



± 



, for 



A = B 
A^B 



(135) 



3 Note that the terms with (5-function are understood in the following way S(x—y) [9(y — x) + 8(x — y)] — 8{x—y). 
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Figure 3: This diagram generates the iV/-dependent part of the term /3q V(x, y) in the two-loop 
ER-BL kernel. 

with 

A GQ F V = 2C F [xlny -x lnx] , A QG F V = 4T F N f -^-_ [xlny - x lnx) , (136) 
A GG F V = 2C A - 2 {y-x). (137) 

y 

In all other cases A AB F(x,y) = 0. 

A consistency check for the /?o-dependent part of the mixed GQ and QG channels was given in 
Ref. |47j by an explicit calculation of the simplest two-loop diagrams with fermion bubble insertion. 
Here we mention that a similar analysis can be easily done for the pure gluonic sector. First we 
should mention that there is no net iV/-dependence of the off-diagonal kernel GG V ND( ' 1 \x, y). 
However, one can argue that among other graphs the diagram in Fig. |3| produces the iV/-part of 
the complete QCD /^-function in Eq. ( |127| ). Since we know that the dotted kernel for the parity 
odd sector, which arises from the logarithmic modification of the leading order kernel, is correct, 
we can just consider the difference between the vector and axial vector channels. Then, neglecting 
in the expression for diagram in Fig. |3| the piece which comes, after subtraction of subdivergencies 
from the 0(e)-effects, (the latter generate only diagonal contributions to the kernel), we get the 
following result for the light-cone position dotted kernel: 

GG JC V - A (y, z) = YlC A yz ln(l -y-z). (138) 

This result is quite obvious, since the only effect of the fermion bubble insertion into the gluon 
line is the shift of the power of the denominator by the amount e. Due to the fact that the weight 
factor associated with this line is (1 — y — z)~ e , the first non-trivial term in the Taylor expansion 
gives us the result fll38|) , where the factor 12C A yz is just the difference between the leading order 
gluonic kernels (j36| , j3^) . By Fourier transformation to the momentum space we easily obtain 

GG V v - A (x,y) = 9(y-x) 

(139) 



A GG F v (x, 



y) + ( GG F v (x,y)- GG F A (x,y) 



v y 6 / 



+ 



x 
JJ 



X 

y j 
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Getting rid of the remaining diagonal piece, we obtain the result displayed in Eqs. ( |135|J136| ). It 
is a trivial task to check that the moments of this kernel is equal to 



( GG 7f 



The last but not the least ingredients for the off-diagonal part of the two-loop evolution kernels 
are the ^-functions, which are given by the formulae 



QQ g(x, y) 
GQ g(x,y) 



-cv 



-C, 



In 1 



0{y - x) 



9{y — x) In 



x 



l- X - 



x 



X 



QG, 



y , 



X 

JJ 



X 

y , 



T g(x, y) = o, 



GG g(x,y) = ^ QQ g(x, y ). 



(140) 



These results were obtained by solving the second order differential equations for the QQ and GG 
channels and by explicit summation of the infinite series of the Gegenbauer polynomials for the 
GQ one. 



8 Solution of the two-loop evolution equations. 

So far we have determined the matrix elements of the anomalous dimension matrix of the conformal 
operators in NLO. In order to solve the renormalization group equation one should determine 
its eigenvectors (the eigenvalues are given by the diagonal elements and coincide with known 
forward anomalous dimensions). This section is devoted to the solution of this problem. It will be 
shown here that they are determined again by the special conformal anomaly matrix introduced 
previously. 

Beyond leading order one has to perform an additional shuffling of operators in order to get the 
eigenstates of the renormalization group equation. Therefore, let us consider the set of operators 

On 

[0;d = £4*(s)[0*], ( 141 ) 

k=0 

where the I?-matrix defines the rotation to the diagonal basis. This new set of operators satisfies 
the renormalization group equation 

^] = -7? PA and l»=[^ G ^y (142) 

with the following solution 

[Ojtifi 2 )] = T exp j-i ^7> s (r))| PM)l (143) 

Since the matrices 7° do not commute with each other, we have introduced the T-ordered expo- 
nential. 



2S 



Substituting Eq. (|141|) into the conformal spin expansion 



<P(x) = J2M^(h'\[0 3J }\h), (144) 

3=0 



for the wave function (j)(x) which is defined as a two-dimensional vector 

'Q<j){x 

and with the conformal waves being Gegenbauer polynomials 



m = Zva ) , (145) 



fa{x) =\ 2 3 2 V , (146) 

we obtain the expansion in terms of the eigenfunctions <pj(x\a s ) of the ER-BL evolution equation 
beyond leading order 

oo oo 

<t>{x) = Y J <l>j{A<Xs){ti\[6 jj ]\h)i with J (x|« s ) = ^0 fc (x) J B fc ,(g). (147) 

j=0 k=j 

As it is seen the eigenfunctions are generalized to non-polynomial functions and corrections to 
the leading order partial waves ( |146| ) are determined entirely by the 5-matrix which satisfies the 
following differential equation: 

/%)|-%) + [7°, B{g)]- + 7 ND £(g) = 0. (148) 

As will be shown in the following two sections, the formalism we have described allows to 
find the corrections to the eigenfunctions analytically, since the I?-matrix is defined in terms of 
the conformal anomaly 7 C and the /3-function. Below we treat two different cases separately: the 
hypothetical conformal limit and the case with running coupling constant. 

8.1 (5{g) = 0. 

As we have observed the problem of finding the corrections to the eigenfunctions is reduced to 
the evaluation of the 5-matrix which satisfies the first order differential equation ( |148| ). In the 
conformal limit of the theory it simplifies to 

[ 7 D , J B(g)]_ + 7 ND J B(g)=0. (149) 

We can construct a recursive perturbative solution of this equation, B = ^=o-^ ■ Since in 
leading order the operators with different conformal spin j do not mix with each other, it provides 
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us with condition B^ = 1. Combining these results give us an inhomogeneous equation for B^ 
with the source 7 ND _B^ -1 ) treated as perturbation 

[7 D , £ w (g)]_ = -j^B^ig), (I > 1). (150) 

For the first non-trivial term B^\ which is the only one we need to the accuracy we are limited 
for, we have the following solution 

*® = -sc^f- (151) 

where on the RHS of Eq. ( |150| ) we insert the expression for 7 ND in terms of the special conformal 
anomaly matrix. Thus, the new two-loop eigenf unctions read 



<j>j{x\a 8 ) = M%) \ hk - 77 -x 7fci \ • (152) 

k=j I a \K,3) ) 

This result can be conveniently written in the form of a convolution, so that we have 

(%) W = {0 '<* -v) + £ (Z Z) s (5) < 153 ' 

Here the matrix-valued function <P(x,y) is expressed in terms of kernels derived in the preceding 
sections 

<2>(x, y) = —{I — V) {S(x, z) ® V® (z, y) + g(x, y)} , (154) 

which thus define the complete set of corrections to the eigenfunctions of NLO evolution equation. 
Here the S'-operator generates the shift in the index of the Gegenbauer polynomials: 



xx] v ^ +p C v j +p {2x-l). (155) 



3 

p=0 



We can also recover the dependence on (3q which is related to the shift of the index of Gegenbauer 
polynomials and not to the effects of a running a s . From the definition of the conformal anomaly 
matrix it follows that this reduces to a mere substitution 

V {0 \x,y) V {0 \x,y) + (j5(x-y). (156) 

Note that the results for the QQ and QG sectors were used by us previously for the evaluation 
of NLO corrections to the amplitudes of the non-forward processes making use of the conformal 
covariant operator product expansion [fy], p^| . 
Let us now address the case with /3(g) ^ 0. 
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8.2 [3(g) ^0. 

In the case of a running coupling constant the solution of the corresponding differential equation 
( |148| ) is more complicated. First, we should specify the boundary condition for B{g). Note that 
in the previous case of a fixed coupling we are naturally led to a scheme with constant B. The 
minimization of the radiative corrections corresponds to the choice B(g ) = 1 that means the 
absence of the former at the reference point fi {go — g'(y u o))- The most important advantage of 
this requirement is that contrary to our previous discussion the initial condition for the solution 
of the RG equation ( |143| ) can be easily determined by composing ordinary Gegenbauer moments 
of the initial wave function since 

PflQk,)] = [OjiMl (157) 
To the required two-loop accuracy we can write the equation for the recursive functions Bjk(g) = 

(3(g)^-B^(g) + [7 D ,B«(g)]_ +7 ND = 0. (158) 
dg 



The solution can be written in the form 

^^Jlm^n^^^Hm^^Y (159) 

where we have introduced ^-ordering along the integration path. 

Recapitulating the results derived here we have found the analytical solution of the two-loop 
ER-BL evolution equation and thus can study explicitly the effects of NLO corrections to the wave 
functions as well as to the non-forward parton distribution functions. 



9 Conclusion. 

To summarize, in this paper we have generalized the formalism, developed previously for the 
evaluation of the two-loop corrections to the non-diagonal part of the singlet exclusive evolution 
kernels in an Abelian theory, to the QCD case. We stress once more that the diagonal entries 
are determined entirely by the known DGLAP splitting functions. An enormous simplification 
of the calculations we have presently performed is due to the fact that the (9(a^)-corrections to 
the evolution kernels are obtained from the evaluation of the one-loop renormalization matrices of 
the conformal operators. To obtain these we have used specific Feynman rules which result from 
the study of the special conformal Ward identities. The latter control the way how the conformal 
symmetry is broken in gauge theories. We would like to point out that our results can be used as 
a guide for a simple convolution-type representation p5| of the full evolution kernels calculated 
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in the direct way. Moreover, numerical simulations previously done in the QQ-channel [[2IJ can 
serve as a consistency check. 

As a by-product of our analysis we gained some insight into the structure of the counterterms 
for the product of two renormalized composite operators, a problem which has not been solved so 
far. Evidently, the complete solution of this task deserves further study in order to put it on the 



same grounds as Joglekar-Lee renormalization theory of composite operators ||57 |. 

We have observed universality of the LO special conformal anomalies in the sense that they 
do not depend on parity and chirality but only on the particles involved. The difference in the 
off-diagonal elements of the NLO anomalous dimension matrices arises as a result of different 
dilatation anomalies, i.e. LO anomalous dimensions of the conformal operators. 

Therefore, at present stage we have all the necessary perturbative input for explicit studies 
of evolution effects for the wave function beyond leading order in the flavour singlet channel. 
However, one should mention that convergence of the partial wave expansion in particular regions 
of the phase space is numerically difficult to achieve. A study of the magnitude of the effects due 
to two-loop evolution for the non-singlet non-forward distribution functions will be addressed by 
us elsewhere [FT . 



Moreover, now it is possible to perform a complete NLO analysis for deeply virtual Compton 
scattering in the leading twist- 2 approximation, since the one-loop coefficient functions were de- 
termined earlier [72], |5T], [73], [711. One can also include the effects of non- leading twist power 



corrections which were already estimated by us using renormalon based techniques [52 
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A Feynman rules for the operator insertions. 

In this appendix we list the Feynman rules for the operator insertions which enter into the special 
conformal Ward identities as well as the rules for the light-cone string operators. The latter when 
sandwiched between appropriate hadronic states define the non-forward distributions or wave 
functions. Namely, these relations read (£ = (p — p')+) 



(h'lWiKuKjlh) 



dx 



^—iK\X — iK,2(C,—X) _ ^—iK2X — %K\{(,—X) 



Q O r (x,(), 



(a.i; 
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(a) 



//, a 
(b) 



Figure 4: Feynman rules for the quark non-local string operator. 
h k 2 



/i, a 



za b 



(a) 




(b) 




Figure 5: Feynman rules for the gluon operator insertions. The grey blob stands either for 
the non-local string operator (<g>), or the z[C^]-vertex (o), or the usual vertices from the QCD 
Lagrangian (•). 



{h'\ G O r {K X ,K 2 )\h) 



— I dx 

2 



e -iKix-iK2(C~x) _j_ e -iK2X-ini((-x) GQ r ( x £\ 



(A.2) 



where the non-local operators are defined in Eqs. ( |7T| , |72D . When ( = 1 and one of the hadronic 
states is replaced by the vacuum the C's coincide with usual singlet wave functions q O(x, ( = 1) = 
x) and g O(x, ( = 1) = G (p(x), which obey the evolution equation ([!]). 
For every gluon and quark line we associate the propagators 



-i)D±{k) 



-iS 



ab 



k 2 + iO 



and 



iS(k) 



k 2 + i0' 



(A.3) 



respectively. The Feynman rules for the vertices are given in the next subsections. The resulting 
expression, composed out of propagators and vertices should be integrated with respect to the 
momentum of every internal line, i.e. multiplied by the factor 



n 



(2vr) M ' 



(A.4) 
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A.l Feynman rules for ®0{k\,K2). 

For the quark string operator we have the following rules 

o = r 



-gt a n.r 



g-iK 2 fc 3+ _ g-iKife 3+ 



(A.5) 
(A.6) 



which are given diagrammatically in Fig. |J Here the upper sign in the square brackets stands for 
r — 7 + , a + j_75 and the lower sign for 7+75. 

A. 2 Feynman rules for g 0(k4, ^2). 

The gluon non-local operator can be decomposed as g O(ki, k 2 ) = g Oa(^i, h 2 ) + g Ona( k i, ^2) + 
k 2 ), where the subscripts designate the origin of the corresponding contributions, namely, 
(N)A: (non-)Abelian part of the field strength tensor; path ordered exponential. In one-loop 
approximation the only Feynman rules we need are the following (crossed circles instead of grey 
blobs in the diagrams of Fig. 0) 



AOf u 



6 ab\ 9*13 ^ f +p . lt (k 1 )f +a . v (k 2 ) e - i ^+-^2+ ± e -^ik2+~i^k 1+ 



(A.7) 



/rtabc 
NA^ pup 



igf 



abc J 9c0 



iK2k 



1 + 



+ n p 



-lKik2+—lK2k3 + 



(At 



e - iKl k 1+ -i K2 k2+ ± ^ K2 )} ; 



*o«% = igf 



abc I 9aP 



(A.9) 



x yi p f +a . u (k 2 )f +/3 . p (k 3 ) 

+ n u f+ a]P (h)f+f3-A k i) 
+ n p f +a Ah)f+pAk2) 



( ,—iKik 3+ -lK2k2+ -p g-«Klfc2+-iK2fc3+ 



^-iK 1 k 1+ -iK2k :j+ __ g— i/ti& 3+ — i/t 2 fel+ 



g— iK\k2+— iK2k\ + _ g— ifti/ci + — iK2^2+ 



fcl + 

g-iK 2 fc2+ _ g-«Klfc2+ 

g-iK 2 fc3+ _ g-«KiA:3 + 
fou. 



Here and below f a /3; P ( k )e p = k a ep — k^e a is the Abelian part of the gluon field strength tensor. 
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A. 3 Feynman rules for z[C^]. 

For the operator insertion i[0^\ which appears in the special conformal Ward identity we have 
(empty blobs in the Feynman rules) 

2 v; b u = -id^u^f^k^^^id^ (ki + fc 2 ), (A.io) 

sV;% = -2g/ afcc {(^i - k 2 ) P 9,u + (h - k 3 ),g vp + (k 3 - ki) u g m } (2vr) 4 2 (h + k 2 + k 3 ^.U) 

4 V;S = 2tg 2 {f abe f cde {g, p g va - g pa g up ) + f ace f bde {g, v9pa - g pa g up ) (A. 12) 

+f ade f cbe {g w g va - g pu g pa )}(27r) 4 2id_5^ (h + k 2 + k 3 + k 4 ) . 

It is instructive to compare them with the usual Feynman rules for the triple and quartic interaction 
vertices (full points in the graphs): 

= gf abc {(h - fc 2 W + ( k 2 ~ **W + (h - hU P ,} (2vr) 4 5( 4 ) (h + k 2 + k 3 ), (A.13) 
= -*g 2 {f abe f cde (9, P 9u, - g^9u P ) + f ace f bde (9,u9 P . - g^gup) (A.14) 

+ f a d e f cbe _ } (2vr )^(4) + k , + k, + fc 4 ) . 

A. 4 Feynman rules for / cfe SB i^ [ G Oji] • 

Finally, in the considered approximation we need only the part of / d d xx^B^(x) SB t^ [ g O] which 
survives for g = 

/ d^xx-B^x)-^^ [ g O{k u K 2 )]\ g=o = ( Kl + K 2 )F^(K 2 n) j/^ j F^imn) (A.15) 

+ (5> 2 n) - (^n)) { 1 i^(«m) - F^{k 2 u) j \ (S«(«m) - <S« ( Kl n)). 

Thus the Feynman rules for this operator are given by («i + k^aO^^i, k 2 ) and the following 
vertex for the second and third term 

iS ab I 9«p 1 _ n * pnii)ha . v{k2 ) + _ W^i)} [e-- lfel +-- 2fe2 + ± e -™ifa+-*"*i + " . 

(A.16) 
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B Renormalization of string operators. 



In this appendix we intend to review the approach for construction of evolution equations for 
non-local light-cone operators developed in Ref. p3| . Although a lot of results have been obtained 
so far using the latter, we did not find in the literature any comprehensive review which allows a 
pedestrian to make these calculations, especially, in the case of non-forward kinematics. To our 



opinion this machinery is more general than the background field formalism ||55|| , since the former 
is not tied to any particular gauge and thus can be used for studying the form of gauge-variant 
counterterms. Next we present the details of the calculations which were cited in the main textQ. 

The momentum integrals which have to be evaluated are of the following form after joining 
the propagators via Feynman parametrization 

f d d q i^^ PjkuX^q) 

J {2n) d [q 2 -L] n ' 1 ' 

where P(ki,Xi\q) is a polynomial function in q, ki and the Feynman parameters Xj. The diver- 
gencies can be evaluated by expanding the exponential factor in the integrand: e - tc i+( K 2~ K i) = 
1 — iq + (K,2 ~ fti) + • • •• Since the denominator depends on q 2 we can average with respect to 
possible rotations of q. Due to the light-like character of the vector n only the first few terms 
survive after integration (we maximally need to expand up to q+). To reduce the result to the 
conventional form, we have to remove the terms proportional to (n 2 — K\) m integrating by parts 
with respect to the Feynman parameters which play now the role of fractions in the light-cone 
position formalism. Typical integration for a test function r(y, z) looks like 

i(«2 - ki)*i+ J * r(y, z) = J* |%)r(0, z) + ^ ^ - 5(1 - y - z)r{y, y)\ , 

and similarly for terms with fc 2+ . In the above equation we have introduced the following shorthand 
notation for the integral 



J= f l dz T dye- ikl+{ y Kl+yK2) - ik2+{lK2+ZKl) . (B.2) 
Jo Jo 



When triple and double integration by parts is required, it turns out that after the first integration 
all dangerous terms proportional to 5(1 — y — z) that might cause a problem vanish identically in 
all cases. 

For the diagrams which originate from the expansion of the phase factor there appear in the 
integrand terms of the type 

E(q, k u y) = r- [l - e -^-^)(yk 1 + q )+] ^ ( R3 ) 



(q + yk 



10 In fact, we have used two different formalisms for our calculations: the one mentioned above and another 
technique developed in Ref. J76|. Computer realizations were also implemented with two different languages: 
FORM and MATHEMATICA, respectively. 
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(a) (b) (c) (d) 

Figure 6: One-loop renormalization of the twist- 2 gluon operator. 

with q being the momentum of integration and y one of the Feynman parameters. The best way 
to treat them is to factorize the g-dependence using the following trivial identity 

1 q+ kl l _ e -i(K2-«i)j/fci+ 

E(q, k u y) = -e~ 9 - . (B.4) 

y h+ 

Expanding as previously the exponential in power series of g™, we can easily perform the final 
momentum integration since only a limited number of lowest order terms contribute due to the 
light-like character of the vector n. 

However, this is not the end of the story since some contributions from the graphs possess 
apart form desired structure k\ + k 2+ also kf + . Therefore, it is necessary to get rid of them by 
reduction kf + k\ + k 2+ . The general structure of these terms in the sum of all diagrams is of the 
form (e.g. the coefficient in front of k\ + ) 

n(y) + r 2 (y)5(z) + r 3 (y)5(l -y-z). (B.5) 

Thus, integration by parts reduces this expression to the desired form up to terms proportional 
to 5(1 — x — y). We get 

k\ + J * n(y) = k 2 1+ J* {[niyjz + T 2 {y)\ 5(1 - y - z) - r 2 (y)5(z)} 
+h+k 2+ J * {[n(y)z + T 2 (y)\ 5(1 - y - z) 

- [ Tl (0)* + r 2 (0)] 5(y) - [r[(y)z + r' 2 (y)]}. (B.6) 

Note also that the resulting function 

k\ + [ Tl (y)z + r 2 (y) + T 3 (y)\ 5(1 - y - z) (B.7) 

can be discard safely because of the Bose symmetry properties of the string operators it is convo- 
luted with. Similar equations hold for k 2+ . 
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Since our result for the GG /C^-kernel is expressed in terms of the leading order gluonic evolu- 
tion kernel we find it useful to give diagram-by-diagram contributions. Namely, straightforward 
calculations of the graphs displayed^] in Fig. || give for parity-odd operators 



GG, 



K {a) (y, Z ) = l-y-z + 4yz+-(l- y)5(z) + -(1 - z)8(y) - (1 - £h5(y)5(z) 



1] 



GG] 



'JC 



(b+c\NA) 



+ h±[2y(l-y)-l} + ^±[2z(l-z 
= -(5 + 0\s(y)5(z) 

+ - y)(2 - y)8(z) + -z)(2- z)5(y), 

A K 2+ A K\+ 



GG] 



'JC(b+ c \$)(y,z 

(d) 



ls(y)6(z) + Uy-3 + 2 



5{z) 



1 



2-3 + 2 



G % d) (y,z) = -yzd(l-y-z) 



' " 1+ -yz5(l - y-z)- \^ ± yz5{\ - y-z), 



(B.8) 

(B.9) 
6(y), (B.10) 

(B.H) 



2/c 2+ "'"" v 9 "' 2k u 

where the symmetry factors are included into the kernels. Making the reduction with Eq. (|B.6| ). 
where T\{y) = 2y(l — y) — l and T2(y) = |(1 —y){2 — y) (and similarly for the y — ► z contribution), 
and summing the resulting contributions with renormalization constant of the gluonic fields 



GGi 



'!C R (y,z) = -^\05(y)6(z), 



(B.12) 



we come to the well known kernel (|86|). 

The contributions to the parity even kernel /C^ from diagrams in Fig. |2] read as follows 



GG, 



TP" 

GG K 



(c 



GGt — 



7C 



GG^- 



(/+</ 



GG^- 



GG: 



(i+fe 
GGjr 



-4i(«i + ^ 2 ) GG /C (a) (2/ ;2 :|e) + W (a+ 6) (B.13) 
-2i(K! + k 2 ) gg /C (o) (</,z|£ = 0) + W (c) (y,z), (B.14) 
2i(«i + «2) GG /C(6+c)(2/ ) *|0 + W (d+e) (y, z), (B.15) 
-2i(«! + k 2 ) { G % b+c) (y,z\0 + G % b+c) (y,z\C = 0)} + W {f+9) (y, z)(B.lQ) 
4i(«! + «: 2 ) GG /C (a) (i/,z|0 + W (/l+i) (y,z), (B.17) 
-4i(«x + « 2 ) GG /C (d) (y, z) + W (j+l) (y, z), (B.18) 

(B.19) 



) = 2i(«! + /€ 2 ) GG /C (d) (i/, z) + WWy, z). 



From these formulae it is clearly seen how the cancellation of gauge dependence in the kernels in 

front of the factor (ki+k 2 ) occurs. To be sure in the correctness of the calculation of the remaining 

part W(y, z), we have kept gauge fixing parameter to be arbitrary. We have observed that the 

11 Of course, there exists a diagram analogous to (d) with composite operator instead of the four-gluon vertex, 
but it is identically zero in dimensional rcgularization since there is no mass parameter entering the loop integral. 
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cancelation of the gauge dependence occurs again between diagrams with the same structure, 
similar to the situation for the kernels which accompany the («i + K 2 )-factor. Namely, we have 



W( a +b+c+h+i)(y,z) 



= —{25(y)-(l-y)8(z)-2yz5(l-y-z)} 
+ - j L{25(z)-(l-z)8(y)-2yz6(l-y-z)} 



W {d+e]NA) (y,z) = ^ jl - 4y + 2y 2 ~(1 +£)<%)} 

+ ^-{i-4« + 2z 2 -i(i + 0*(«)}*(y) ) 



(B.20) 



(B.21) 



W( d +e\t>){y,z) 



+ 



2^-3-2 
2y - 3 - 2 



W(f+g\NA)(V,z) 

W( f+g \<i>)(y,z) 



1 

h+ 
1 

k l+ 
1 

- 1 - 1 1 - 2z + 6 

^2+ 1 



+ 2 



+ 2tf(s) <%) 



+ 2%) 



(B.22) 



3y - 2y l + \b{y) S(z) + — l3z-2z 2 + ±6(z) %),(B.23) 



+ 



1 

fen 



1 - 2y + 6 



fci 



(B.24) 



W(j+k)(y,z) 

yv { i)(y,z) 



^- {Ayz -\}5{l-y-z) + ^- {Ayz -l}5(l-y-z), (B.25) 
fc 2 + fci+ 

{1 - 2y*} 5(1 - y - z) + ^- {1 - 2y*} 5(1 - y - z). (B.26) 



To simplify the result given by W( a +b+c+h+i)(y, z ), we have used Eq. (P-6Q and the following 
identity 

J * (h+ + h+)6(l ~ y -z)=J* {k 1+ 5(z) + k 2+ 5(y)} . 



(B.27) 



Finally, the kernel which we are interested in is given by 



GG) 



1 



K w (y,z) = W(y,z) + 2i— + - 



Ca 



1+ ^2+ 
1 ■ 
72 



- 2 



8{z) 



(B.28) 



where we restored the colour factors omitted previously. This kernel arises as a factor convoluted 
with J . Of course, it is easy to integrate by parts to get rid of momentum factors ki + and acquire 
instead i(K 2 — Adding then the contribution with k 2 *-> «i gives the resulting kernel in front 
of the Bose symmetrical gluon factor ( A.7| ) . However, we will skip this last trivial step since we 
can obtain the desired momentum fraction ^-function from equations similar to ( B.28|) with the 
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substitution^ = y, k 2 + = y and Fourier transformation. 

Finally, let us mention that the results of the calculations for the parity odd gluonic sector 
are the same as given above in Eq. ( |B.28|) . The differences in the intermediate steps arise from 
the fact that the contributions of the type yz5(l — y — z) are absent by simple reason of Bose 
symmetry of the light-ray operators which are convoluted with this kernel. 



C Properties of hypergeometric functions. 



In this appendix we list all formulae which are necessary for the evaluation of the Gegenbauer 
moments of the evolution kernels. 

The generalized hypergeometric function is defined as infinite series [[75fl 

Ot\j ■ ■ ■ , Oi m 

" ~&(iM/?i)<- ■■(&)/ 



F 



X 



(C.l) 



Pi,...,P„ 

where we have introduced the Pochhammer symbol (a)f = T(a + £)/T(a). The integral represen- 
tation of the hypergeometric function is (where iF (a\xy) = (1 — xy)~ a ) 

a m+ i 



X 



Pi, ■ ■ ■ , Pn+1 

r(/Wi 



(C.2) 



— f 1 dyy a ^- l y am+1 - Pn+1 - l m F n ( ' ' ' 
-i) Jo \Pi,--- 



OL\y • • • , 0( m 

,Pn 



xy 



T(a m+ i)T(p n+1 - a r 

When one of the upper indices otj = —j, where j G IN, the infinite series ( |C.1|) degenerates 
into a finite polynomial of order j. Therefore, the relation between the Gegenbauer polynomials 
and the hypergeometric function reads 



X 



(C.3) 



In the calculation of the integrals we use the following representation of the Gegenbauer poly- 
nomials (w(x \v) 



(xxY 2 



C v A2x 



•,2v 



r(i)r(z/ + J + i) 1 } dx> 



'xxV +U -2. 



(C.4) 



Integrating j-times by parts we end up with j derivatives acting on the hypergeometric function. 
This differentiation can be done with the help of the formula 



-x 



fe 771 / ' ' ' - J 



X 



dxi \ Pt, . . . ,p n 

(oii)j-fc ■ ■ ■ {a m )j-kj ] - 

(Pl)j-k ■ ■ ■ (Pn)j-k(j ~ k)\ 



"l + j - k, 
Pl+j-k,.. 



.,a m +j-k,j + l 
P n +j-k,j-k + l 



(C.5) 



for j — k > 0. 



12 Do not mix this y variable (which is the momentum fraction in the ER-BL kernel) with y in the light cone 
position formalism which is integrated out with Fourier transformation. 
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To simplify the result of integration it is instructive to use a relation between the functions of 
(m, n) and {m — 1, n — 1) orders 



a 1 ,...,a r 



Pi) ■ ■ ■ , Pn-lj a m + 1 



.r 



Q! m / Ot\, • • • j 0£ m —i 

l + a m -f3 n -i m n \ Pi, . . . ,(3 n -i 



X 



(C.6) 



Pn- 



1 



F 

m ± n 



Oil, ■ ■ ■ ; C^m 
P\i- ■ -i Pn-1 ~ 1) a m + 1 



.r 



1 + «m - Ai-l ' 

The number of consequent iterations is governed by the condition of fulfilling of the relation 
1 + YU a i — Pj- Once being satisfied it allows to apply Saalschiitz's theorem []75|] 

\ _ (j-p^l + a-j^ 



3 r 2 



+G!,P 

7, l + a + /3-7 



(iW + ^ + P-iV 



(C.7) 



From this relation follows immediately 



3-^2 

For 2 F\ we have 
2F1 



-j,j + a,p 
a - 1, 1 + P 



P 



a — 1 fa) 



-iy + 



a — P — 1 (a — 

+ 



./ J 



(C.8) 



/3 



1 



1 



.(l + a-/fy 



-j + k,j + k + 3 
2k + 4 



1 



(C.9) 



Since for the regularized functions we get the derivatives of the hypergeometric function 3F2 with 
respect to an index, we need the following expansion 



3 r 2 



-j,j + <x,P 
l + a,P + e 



r(l + q) 

r(i + a+i; 



^■o + er(j) 



0(e 2 ), (CIO) 



which can be derived making use of a fundamental identity for 3 .F 2 -functions [[75|] and the definition 
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